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Many structural materials rely on the formation of chromia, silica, or alumina as a 
protective layer when exposed in high temperature oxidizing environments. Growth of 
these oxides provides a protective diffusion barrier which slows down further oxidation. 
In water vapor, however, each of these oxides can form volatile hydroxide species which 
remove the surface oxide, thus, lowering the protective capability of the oxide scale. 
This paper summarizes the current understanding of volatility of chromia, silica, and 
alumina in water vapor-containing combustion environments. In addition unanswered 
questions in each system are discussed. 

It is well known that chromia scales volatilize in oxidizing environments to form 
Cr 03 (g). In water vapor, the very stable Cr02(0H)2(g) species is formed which has a 
much higher partial pressure than CrO.^g). However, there are at least three sets of 
thermochemical data for this species which differ by as much as two orders of magnitude 
in predicted partial pressures of Cr02(0H)2(g). Differences between these data sets have 
not yet been resolved. 

Within the last ten years, the importance of silica volatility in combustion 
environments has been identified. Silica is known to react with water vapor to form 
Si(OH) 4 (g). Accurate thermodynamic data are available for this species. However, there 
are a number of unanswered questions in this system. First, at high temperatures 
(>1300°C), there is evidence for SiO(OH) 2 (g) formation from silica and water vapor. 
However, available thermodynamic data for this species show large variations. Current 
work is underway to resolve this issue. Second, silica volatility rates in fuel-rich 
combustion environments are more rapid than predicted for Si(OH) 4 (g) formation alone. 
A third unanswered question is why silica volatility rates appear to differ depending on 
whether silica is formed on silicon carbide or silicon nitride in combustion environments. 
One possible explanation is that the silica volatility is the same, but differences in 
material temperature at a constant gas temperature result in apparent differences in 
volatility. Finally, at very high gas velocities silica is absent from the surface of silicon 
nitride, while recession rates are consistent with silica volatility by Si(OH) 4 (g) formation. 
Possible active oxidation-water vapor reaction mechanisms are discussed. 

Alumina volatility in water vapor has also been recently observed, albeit at higher 
temperatures than chromia and even silica volatility. It is believed that the formation of 
Al(OH) 3 (g) is the primary volatilization reaction in combustion environments. The 
identity of this molecule has been proposed based on the observed water vapor pressure 
dependence of its formation. Preliminary results showing mass spectrometric 
identification of this species will be presented, but more definitive results are needed. 

In conclusion, a comparison of the three systems shows that the partial pressures 
of volatile metal hydroxides formed in combustion environments decrease in magnitude 
in the following order: Cr02(OH) 2 (g)>Si(OH)4(g), SiO(OH) 2 (g)>Al(OH) 3 (g). 
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Introduction 

Many structural materials rely on the formation of chromia, silica, or alumina as a protective layer 
when exposed in high temperature oxidizing environments. The presence of these oxide layers 
provides a protective diffusion barrier which slows down further oxidation. In atmospheres 
containing water vapor, however, reactions to form volatile hydroxide species occur which remove 
the surface oxide, thus, lowering the protective capability of the oxide scale. This paper 
summarizes the current understanding of volatility of chromia, silica, and alumina in water vapor- 
containing combustion environments. In addition unanswered questions in each system are 
discussed. These systems were reviewed in a previous publication [1]. The current paper 
represents an update on the considerable information learned in the past five years for these 
systems. 

Volatility is defined here as material loss due to the formation of vapor species. Volatility of 
oxides in environments containing water vapor occurs by the following generalized reaction: 


M0 X + nH 2 0(g) + m0 2 (g) = M 0 (x+ n+ 2m) H (2n) (g). ( 1 ) 

For the systems of interest, here, the specific volatilization reactions are: 

>/ 2 Cr 2 0 3 + H 2 0(g) + % 0 2 (g) = Cr0 2 (0H) 2 (g) (2) 

V4 Cr 2 0 3 + 3 / 4 0 2 (g) = Cr0 3 (g) (3) 

Si0 2 + 2H 2 0(g) = Si(0H) 4 (g) (4) 

Si0 2 + H 2 0(g) = SiO(OH) 2 (g) (5) 

Si0 2 + l / 2 H 2 0(g) = l / 4 0 2 (g) + SiO(OH)(g) (6) 

l / 2 A1 2 0 3 + 3/2 H 2 0(g) = Al(OH) 3 (g) (7) 


For systems in which more than one volatilization reaction occurs, the reaction listed first is of 
primary interest for typical temperatures and water vapor partial pressures of industrial processes 
such as cooling, processing, and combustion applications. These reactions will be discussed in 
more detail in a following section. 


Kinetic Models 

This paper is concerned with the rate at which these volatilization reactions occur and the resulting 
consumption rate of underlying alloy or substrate. As a result, kinetic models which predict these 
rates are described next. 

Tedmon oxidation/volatilization kinetics. At short times or for moderate oxide volatility, the 
oxide thickness is a function of both the parabolic oxidation rate, k’ p , and the linear volatility rate, 



k’l (in units of length and time). Tedmon [2] has described oxidation/volatilization kinetics for 
chromia formers using the following relationship: 
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where x is the oxide thickness and t is exposure time. So that after integration the relationship 
between exposure time and oxide thickness can be expressed as: 
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At long times or high volatility rates, a limiting oxide thickness, xl, is developed which is given by: 


X L 



( 10 ) 


Similar expressions can be developed for weight change and material recession. Thus the oxide 
thickness can be calculated if the rate constants are known. 

Measured Oxidation Rate Constants. The oxidation rate constants are fairly well known for 
chromia, silica, and alumina formers. Chromia formation rates have been summarized by Hindam 
and Whittle [3], There is large scatter in the oxidation rates, with those alloys containing reactive 
elements showing lower oxidation rates. The oxidation rate for chromium [4] and for chromia- 
forming alloys were found to be independent of oxygen partial pressure [3,5]. It has not been 
demonstrated that the oxidation rate is independent of water vapor partial pressure, although it is 
assumed to be true for the purposes of this paper. 

The oxidation rates of Si [6], SiC [7,8], and SriN 4 [9] have been measured in water vapor. Silica 
formation rates are known to depend on the pressure of the oxidant as well as the type of oxidant. 
It is generally agreed that amorphous silica growth has a pressure dependence described by a power 
law exponent of one, or close to one, in both dry oxygen and water vapor [6,8], The oxidation rate 
of silica formers is higher in water vapor than in dry oxygen since water vapor has a higher 
solubility in amorphous silica compared to oxygen [6], 

Oxidation rates for alumina-formers have also been summarized by Hindam and Whittle [3]. 
The measured oxidation rate of NiAl [10] is considered here as representative for pure alumina 
formers. Alumina formation rates from an alloy have been shown to be independent of oxygen 
partial pressure [11]. Oxidation of alumina-forming alloys in the presence of water vapor does not 
change the kinetics unless the oxidation is cyclic in nature [12,13], 

Calculated Volatilization Rates. Volatilization rates can be calculated once the gas flow 
conditions over the oxide are defined. In many combustion environments, volatility is limited by 
transport of the volatile species away from the oxide surface through a gas boundary layer. To 
calculate the volatility for the boundaiy layer limited process, the equilibrium partial pressure of the 
volatile species must be known, or thermodynamic data must be available to calculate this pressure. 
Sources and limitations of the available thermodynamic data for the systems of interest in this paper 
are reported in the next section. Assuming volatility is limited by transport of the volatile species 
through a laminar boundary layer over an oxide of flat plate geometry, the following expression can 
be used to calculate the volatility rate [14,15]: 



( 11 ) 
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where ki is the flux of volatile species in g/cm 2 h, Re is the dimensionless Reynold’s number, Sc is 
the dimensionless Schmidt number, p is the density of the volatile species in the laminar boundary 
layer, D is the interdiffiision coefficient of the volatile species in the boundary layer and L is the 
characteristic length of the flat plate. The interdiffusion coefficient is calculated using the 
Chapman Enskog equation [15] and the collision diameter and collision integral for the volatile 
species. These terms are generally not known and must be estimated from molecules of similar 
structure. The rate, ki, is given in terms of weight loss of oxide in Eq. 1 1 and is related to the rate, 
k’i, given in terms of oxide recession in Eqs. 8 through 10 by the oxide density. Expanding the 
terms for Re, Sc, and p results in the following expression: 
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where p’ is the gas density of the boundary layer, v is the linear gas velocity, rj is the gas viscosity, 
M is the molecular weight of the volatile species, R is the gas constant, T is absolute temperature, P 
is partial pressure of the noted species, n and m are the pressure dependence exponents for 
formation of one molecule of volatile hydroxide as in Eq. 1, AG is the free energy of reaction for 
formation of one molecule of hydroxide from the oxide and water vapor as in Eqs. 2 through 7, and 
all other terms have been previously defined. The activity of the solid oxide is assumed equal to 
one in Eq. 12. Simplifying Eq. 12 in terms of pressure (total pressure, P, as well as partial 
pressures of water vapor and oxygen), gas velocity and temperature results in the following 
relationship: 
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Here, the temperature dependencies of p’, q, and D are neglected, and a simple Arrhenius 
temperature dependence is assumed based on the free energy of reaction to form the volatile 
hydroxide. At high velocities, when flow is turbulent (Re>3xl0 5 ), the Reynolds number will have 
an exponent of 0.8 rather than 0.5 [14] and Eq. 12 should instead be simplified to: 


kj ocv 


p n p m 
0.8 H 2 0 0 2 


,0.2 


exp - 


AG 

RT 


(14) 


The thermodynamic data needed for the calculation of volatility rates are described in the next 
section. 


Thermodynamic Data for Volatilization Reactions of Common Oxides 

The chromia-water vapor system. Three sets of thermodynamic data exist for the Cr02(OH)2(g) 
species. Ebbinghaus [16] evaluated the data available in the literature up to 1993. Gindorf et al. 
[17] have measured equilibrium partial pressures of CrC> 2 (OH) 2 (g) under typical solid oxide fuel 



cell cathode conditions using a transpiration technique. Finally, the IVTAN database [18] also 
contains data for this species, although the method of determination of these data is unknown to this 
author. Fig. 1 shows the disparity in calculated partial pressures of Cr 02 ( 0 H )2 based on the three 
data sets. Clearly, the differences in these data must be resolved before accurate predictions of 
chromia volatility can be made. C 1 O 3 is also an important volatile species at low water contents 
and high temperatures. Data for this species from a number of sources are reliable [16,19,20]. 
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Fig. 1. Variation in calculated and measured pressures of Cr 02 ( 0 H )2 at 0.02 bar H 2 O and 0.21 bar 
O 2 from three sources of thermodynamic data. Data for CrC >3 are also shown. After Ref. [17]. 

The silica-water vapor system. Four sets of thermodynamic data are known for the Si(OH) 4 (g) 
species. In this case, the results are in good agreement. Krikorian [21] has estimated 
thermodynamic data based on partition functions and comparison to molecules of similar structure. 
Hashimoto [22] has measured equilibrium partial pressures of Si(OH) 4 (g) using a transpiration 
technique. Allendorf [23] has calculated thermodynamic data for this species using ab initio 
techniques. Finally, Copland et al. [24] have also measured equilibrium partial pressures of 
Si(OH) 4 (g) using a transpiration technique. Experimentally measured volatility rates [25,26] are 
consistent with those rates calculated with Eqs. 13 and 14 and thermodynamic data for Si(OH) 4 (g). 
Other volatile Si-O-H species, such as SiO(OH )2 or SiO(OH) may be important at temperatures of 
1400°C and higher. Thermodynamic data for these species are not considered reliable at this time. 

The alumina-water vapor system. Three sets of data are available for the Al(OH) 3 (g) species. 
Gurvich et al. [20] have estimated thermodynamic data for this species based on partition functions 
and the structure of similar molecules. Hashimoto [22] has measured equilibrium partial pressures 
of Al(OH) 3 (g) from a mixture of calcium aluminate phases in water vapor using a transpiration 
technique. Finally, Allendorf [27] has calculated thermodynamic data for this species using ab 
initio techniques. All three sets of data are in fairly good agreement. Experimentally measured 
volatility rates [26,28] are consistent with those rates calculated with Eqs. 13 and 14 and 
thermodynamic data for Al(OH) 3 (g). 

Comparison of chromia, silica, and alumina volatility. The relative rates of volatility for these 
oxides in water vapor can be assessed by comparing the equilibrium partial pressures of volatile 
species formed for a given set of conditions. Partial pressures for the primary volatile species in 
each oxide system have been calculated at 1 bar total pressure with equal amounts of oxygen and 
water vapor over the temperature range of 700 to 1500°C. The comparison is shown in Fig. 2. 
Both the Ebbinghaus data [16] and the IVTAN [18] data for Cr 02 ( 0 H) 2 (g) are plotted to span the 


possible range of partial pressures for this species. CrOa(g) is neglected on this plot although it 
would become important under these conditions at temperatures greater than about 1100°C if the 
IVTAN values for Cr02(0H)2 are correct. Gurvich data [20] are used for Al(OH) 3 (g) while 
Hashimoto’s data [22] are used for Si(OH) 4 . Other Si-O-H species may become important at 
temperatures above 1400°C. 
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Fig. 2. Partial pressures of primary volatile species calculated for each oxide at one atm total 
pressure: 0.5 atm O 2 , 0.5 atm H 2 O. 


Mapping Recession in a Combustion Environment as a Function of Total Pressure and 
Temperature 

Recession from steady state volatility can be determined as a function of pressure, temperature, and 
gas velocity as described by Eqs. 13 and 14. The velocity dependence of volatilization will be the 
same for all oxide systems in a given application since it depends only on the flow characteristics of 
the gaseous environment. In combustion environments, the oxygen partial pressure and water 
vapor partial pressure will scale with the total pressure. For fuel-lean hydrocarbon combustion 
(equivalence ratio=0.6), the oxygen partial pressure and water vapor partial will be about 8% and 
9% respectively of the total pressure [29]. The pressure dependence can thus be simplified by 
replacing Ph 2 o by 0.09 P to tai and P02 by 0.08 P to tai leaving one pressure term in Eqs. 13 and 14. 
Variation of recession limited lifetime in combustion environments as a function of pressure and 
temperature can be visualized with the aid of a recession map as shown in Fig. 3 for Cr02(0H)2 
volatilization from chromia. Again, both the data of Ebbinghaus [16] and IVTAN [18] were used 
to calculate steady state volatility rates at a gas velocity of 50 m/s as a function of temperature and 
total pressure in a combustion environment. The pressure-temperature conditions are plotted which 
show a recession of 250 pm in 1,000 or 100 hours. The Ebbinghaus data are plotted for lOOh, since 
lOOOh lifetimes could only be met for total pressures less than 1 atm. 

Fig. 4 compares oxide growth rates to oxide volatilization rates at 10 atm total pressure with a 
gas velocity of 50 m/sec. The oxidation curve is based on the expression for k p from Hindham and 
Whittle [3] for oxidation limited by transport of interstitial chromium ions. This rate is high 
relative to oxidation of chromia formers with reactive elements. Overall, the recession of chromia 
by volatility compared to recession from oxidation shows the relative importance of chromia 
volatility in combustion environments. Note that this comparison plot again neglects the effects of 
Cr 03 (g) formation which would be significant if the IVTAN data for Cr02(0H)2(g) are correct. 


Similarly, recession maps are shown for silica and alumina in combustion environments in Figs. 
5 and 6 respectively. One obvious result of these maps is that the alumina-water vapor system is 
much less sensitive to recession by volatility compared to the other oxide systems. Only at high- 
temperatures is alumina volatility an issue. Thus, alumina volatility has not been a problem for 
alumina-forming superalloys since their use temperature is below 1200°C. Water- vapor induced 
volatility will be more of a concern at the higher use-temperatures proposed for alumina-based 
oxide/oxide composites and alumina-containing environmental barrier coatings. 



Fig. 3. Recession map for chromia volatility 
at a gas velocity of 50 m/sec. Pressures and 
temperatures below the curves represent 
conditions with acceptable recession rates. 



Fig. 5. Recession map for silica volatility at a 
gas velocity of 50 m/sec. Pressures and 
temperatures below the curves represent 
conditions with acceptable recession rates. 
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Fig. 4. Comparison of times to grow or 
volatilize a chromia scale of 250 pm as a 
function of temperature. Volatility calculated 
for P = 10 atm and v = 50 m/sec. 



Fig. 6. Recession map for alumina volatility. 
Pressures and temperatures below the curves 
represent conditions with acceptable 
recession rates. 


Unanswered questions for chromia, silica and alumina volatility in high-temperature water 
vapor 

The chromia-water vapor system. Clearly, the biggest question in this system is the stability of 
the Cr02(0H)2(g) species. The discrepancy between the Ebbinghaus [16], Gindorf [17], and 
IVTAN [18] data must be resolved before accurate chromia volatility predictions can be made. 
Weight change measurements of 304L stainless due to oxidation/volatilization in water vapor- 
containing environments by Asteman et al. [30] indicate that the Ebbinghaus data over-predict the 
volatility of chromia in water vapor. 


The silica-water vapor system. This system has been studied more thoroughly at NASA-Glenn, 
thus the higher familiarity with this system has raised more questions. First, and most importantly, 
at temperatures greater than 1400°C, there is evidence that an Si-O-H species, in addition to 
Si(OH) 4 , contributes to silica volatility. Preliminary evidence suggests that SiO(OH )2 may explain 
these results [31]. The transpiration study by Hashimoto [22], however, showed no evidence of any 
species other than Si(OH )4 to temperatures as high as 1500°C. 

Similarly, silica volatility in fuel-rich conditions can not all be explained by Si(OH) 4 (g) 
formation. The additional volatility could be explained by an Si-O-H species with a higher water 
vapor partial pressure dependence titan Si(OH) 4 (g) [32,33]. 

Differences in silica volatility were found for SiC and Si 3 N 4 in burner rig exposures [34] even 
after correction for volatilization of CO or CO 2 and N 2 . One possible explanation is that the 
difference in emissivity of the materials results in different surface temperatures in the same 
environment. 

Finally, Tedmon kinetics and volatilization limited by transport through a gas boundary layer no 
longer accurately describe the volatilization mechanism when gas velocities and the corresponding 
volatilization rates are so high that no silica is found on the surface of a silica-forming material. 
Bare Si 3 N 4 surfaces have been observed in turbine vane testing [35,36]. In this case, the 
volatilization and recession could be limited by oxide growth, or a completely different mechanism 
similar to active oxidation could be operative. A mechanism can be envisioned in which SiC or 
Si 3 N 4 volatilizes directly to Si(OH )4 without the formation of a solid silica intermediate step. The 
temperature dependence of such a process is expected to vary from that of the 
oxidation/volatilization process. 

The alumina-water vapor system. Fairly good agreement exists between estimated data for 
Al(OH) 3 (g) formation, transpiration results for a mixture of calcium aluminate phases, weight loss 
of sapphire coupons, and volatility observed in a high pressure burner rig. The major shortcoming 
of the data for this system is that the existence of the Al(OH) 3 (g) species has not been confirmed 
experimentally by any technique, but only inferred from the pressure dependence of volatility. 
Preliminary mass spectrometry results have been obtained [37] for Al(OH) 2 + at m/e=61 amu, a 
possible fragment of Al(OH) 3 , but more definitive evidence is required to confirm the identity of 
Al(OH) 3 (g). 


Summary 

The volatility of chromia, silica, and alumina in high-temperature water vapor has been discussed. 
Kinetic equations which describe both the oxidation and volatilization processes have been 
reviewed. Recession maps have been presented for each oxide-water vapor system which show 
pressure-temperature use ranges for a given recession limit. Recession rates have been compared to 
oxidation rates. Topics related to oxide volatility which require further understanding have been 
summarized. Finally, it has been shown that the severity of volatility increases in order for the 
alumina-, silica-, and chromia-water vapor systems. 
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